Casimir Effect for Massless Fermions in One Dimension: A Force Operator Approach 
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We calculate the Casimir interaction between two short range scatterers embedded in a back- 
ground of one dimensional massless Dirac fermions using a force operator approach. We obtain the 
force between two finite width square barriers, and take the limit of zero width and infinite potential 
strength to study the Casimir force mediated by the fermions. For the case of identical scatterers 
we recover the conventional attractive one dimensional Casimir force. For the general problem with 
inequivalent scatterers we find that the magnitude and sign of this force depend on the relative 
spinor polarizations of the two scattering potentials which can be tuned to give an attractive, a 
repulsive, or a compensated null Casimir interaction. 
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Boundaries modify the spectrum of zero point fluctua- 
tions of a quantum field, resulting in fluctuation-induced 
forces and pressures on the boundaries that are known 
generally as Casimir effects [ij. When sharp bound- 
aries conditions are used to model the Casimir effect, 
they yield perfect reflection of the incident propagat- 
ing quantum field at all energies [ij. However, in many 
physical applications this hard-wall limit is not appro- 
priate; of special interest in the present work are interac- 
tions between localized scatterers in one dimension that 
have energy-dependent scattering properties controlled 
by the strength, range and shape of the potential. Along 
this line, previous work has recognized that the finite 
reflectance of partially transmitting mirrors provides a 
natural high energy regularization scheme for comput- 
ing the effect of sharp reflecting boundaries on the zero 
point energy of the electromagnetic fleld 0, Q • In more 
recent work, Sundberg and Jaffe approached the problem 
of computing the effect of conflning boundary conditions 
on a degenerate gas of fermions in one dimension as the 
limiting behavior for rectangular barriers of flnite width 
and height. Interestingly, they encounter a divergence 
of the Casimir energy in the zero width limit (a sharp 
boundary) even for flnite potential strength 0]. 

In this Rapid Communication we address the problem 
of Casimir interactions between scatterers mediated by a 
one-dimensional Fermi gas. The fermions in our calcula- 
tion are massless Dirac fermions appropriate to describe, 
for example, the (single-valley) electronic spectrum of a 
metallic carbon nanotube. We employ the Hellmann- 
Feynman theorem to calculate the force, rather than en- 
ergy, of interaction between two scatterers as a function 
of their separation d. This approach renders our calcu- 
lation free from ultraviolet divergences even for the lim- 
iting case of sharp scatterers. We demonstrate that for 
the case of identical scatterers, this formalism recovers 
the well known attractive l/dP' Casimir force in one di- 
mension. Furthermore, we flnd that for Dirac fermions 
the internal structure of the matrix-valued scattering po- 



tential admits a long range Casimir interaction which can 
also be repulsive or even compensated. This provides a 
physical situation where the Casimir interaction is con- 
tinuously tunable from attractive to repulsive by varia- 
tion of an internal control parameter, realizing the known 
bounds for the one dimensional Casimir interaction as 
two limiting cases. The results may be relevant for in- 
direct interactions between defects and adsorbed species 
on carbon nanotubes. 

The fermions in our model are massless one- 
dimensional Dirac fermions described by the Hamiltonian 

iy~iaA + V{x)-E)^k{x)={), (1) 

where we set h ^ c = 1. In graphene and carbon nan- 
otubes the spinor polarizations describe the internal de- 
grees of freedom generated by the two-sublattice struc- 
ture in its primitive cell. When V{x) ~ 0, the eigen- 
states of Tio are plane waves multiplying two-dimensional 
spinors, '^k{x) — ^k^-^^^ /V^- When the chemical po- 
tential is flxed at /i = 0, the fllled Dirac sea has E = —\k\ 
with = (l,Tl)/V2. 

The general form of the potential entering (1) is 
V{x) — Vo{x)i + V{x) ■ a. The part of the poten- 
tial can eliminated by a gauge transformation and a 
scalar potential proportional to the identity matrix pro- 
duces no backscattering in the massless Dirac equation. 
Therefore, we consider potentials for which V lies in the 
2/z-plane. In this paper, we consider the effects of the ori- 
entation of the potential determined by angle 0. Thus, 
a square barrier potential located between points xi and 
X2 is written as 

V{x,<j))^V{(b)0{x~xi)eix2-x), (2) 

where V{(p) = Ye^<y.4>/2^^^-tT.ck/2 ^ ^^^^ q(^^-^ ^ g^^p 
function. 

To study the force on a square well scatterer we use 
the Hellmann-Feynman theorem, (07i(A)/9A) = dE/dX 
Taking the control parameter A = {xi + X2)/2 = x, 
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the ground state average gives the force acting on a rigid 
barrier. For a barrier with sharp walls the expectation 
value becomes 

(*(:^)| ^ = {^{x + a/2)\ V + a/2)) 

-{^ix-a/2)\V\^{x-a/2)), (3) 

where V is the square barrier potential, x is its center and 
a is its width. The total force is the expectation value of 
this force operator, F — —dH/dx, summed over all the 
occupied states; Eq. ^ then gives the difference between 
the pressures exerted on the right and the left sides of 
the barrier. For potentials of general shape a similar 
expression can be developed in terms of an integral over 
the scattering region. 

First, we apply Eq. ^ to calculate the force on an iso- 
lated barrier. The eigenstates are represented as linear 
combinations of right and left moving solutions of Tio'. 
^'(a;) = -^{ak^ks''''' + /3fc$-fce-''=^), where and 
(3k represent the amplitudes of the counterpropagating 
waves in each region. The yz polarized potential defined 
in Eq. ^ gives V{(p)<i>±k = Ve^*'l'<S>^k, so the general 
expression for the expectation values in Eq. ([3]) at some 
position X is 



V, 



{^{x)\V{<l>)\^{x)) = -7^e[afe/?^e^(2fex+« 



(4) 



We use a transfer matrix to obtain the coefficients ak 
and Pk entering Eq. (jl]). The transfer matrix is defined 
so that ^{x2) = T'i'{xi), where xi and X2 are the left 
and right boundaries of a barrier, respectively ; T is cal- 
culated by integrating Eq. [l] 



T ^ Px exp 



dxax[E ~V{x)] 



(5) 



where Px is a spatial ordering operator. For the square 
potential of width a defined in Eq. ([2]) , the transfer ma- 
trix for negative energy states is 



T = cos((7a) 



[X X V) 



sin(ga), (6) 



where V = V^(0, sin(/>, cos i/i) defines a potential in the 
yz-plane, q = — V^, and fc > 0. 

From the transfer matrix we calculate the scattering 
matrix 5, which gives the transmitted (t) and reflected 
(r) amplitudes for wave incident on the barrier from the 
right and from the left. The unitary S-matrix for a single 
square barrier is 



^1 



re 



i{2kxi+4i) 



re 



-i{2kx2+<t>) 



te 



-ika 



(7) 



The transmission and refiection coefficients can then be 
parameterized t = re"' and r — — r^e"', where 



A 



(1/2 cosh^ Aa-/c2)i/2 



, rj — tan 



-1 



k tanh Aa 
A 



(8) 



with A = —iq = \/V'^ — fc2. To obtain the hard-wall 
limit, we fix r = Va^ and take F oo. In this limit, 
|r|2 1 and jip at all energies. 

For a single barrier, the contributions to the force from 
the particles incoming from the right and the left cancel, 
resulting in no net force. A nonzero force arises from the 
multiple reflection of electron waves between two barri- 
ers. An illustration of a scattering process for two square 
potentials with different spinor polarizations 4>i and 02 
separated by distance d is shown in Fig. [T] The con- 
tributions from waves incoming from the right are also 
included in the calculation. 
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FIG. 1: Scattering of massless Dirac fermions (incoming from 
the left) between two square barriers of height V, width a, 
and separation d. The two potentials defined in Eq. ((2)1 have 
a spinor polarization determined by angle 0. The reflection 
and transmission coefficients are labeled in each scattering 
region. 

The S-matrix for the two-barrier system fel in Fig.[T] IS 

S - f ^ ^""'^'^ (Q) 



The total reflection and transmission coefficients shown 
in regions I and III of Fig. [l]are given by 



T = ^— , R = re-^'^i^^+d) (i + L 

l-r^e'" V 1 



^2gi(2fca+iy) 



, (10) 



where v — 2kd -\- dcj) and 5(j) — (1)2 — 4>i. Ti and i?i in 
region II of Fig. [1] are given by 



1 — r^e 



, Ri 



1 — r^e 



2 f^iu 



(11) 



The coefficients for the waves incoming from the left (i?i 
and Ti), and the ones incoming from the right [R'l and 
T[) are related by R{ = i?ie-'('^i+'^2) and T[ = Ti. 

To calculate the force in the two-barrier problem we fix 
the position of the left barrier in Fig. [T] and differentiate 
the Hamiltonian with respect to d. To obtain the total 
force, we sum over the occupied states of the filled Dirac 
sea at fixed chemical potential. We find that the force 
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between two square barriers of finite height and width is 
Jo 2^ 

-RiT*e"'^'"^+'l'^\l + e''')]. (12) 

The first term in the integrand arises from the exterior 
modes pushing the two barriers together. The second 
term accounts for the confined modes in between the bar- 
riers pushing them apart. Since incoming waves are fully 
transmitted at high energies for barriers of finite height 
and width, the integral in Eq. converges even in the 
case of sharp barriers (a ^ 0), with F = Va fixed. Thus, 
the refiection coefficient provides a natural cutoff for the 
computation of the force (though not the energy [3|) even 
in the limit of infinitely high barriers. 

The Casiniir force for hard-wall boundary conditions 
requires the limits of infinite barrier strength F — s- oo 
and zero width a — > 0. This limit enforces a vanishing 
current at the boundaries, the so-called bag boundary 
conditions. Since the force in Eq. (fT2|) is multiplied by 
V, we keep terms to 0{k/V) in the integrand. The first 
term in Eq. ()12|) becomes proportional to fc, thus imply- 
ing a continuous spectrum of modes scattering off the 
barriers from the outside. The second term exhibits res- 
onances that arise from the quantized modes between 
the boundaries. These resonances, similar to ones seen 
in Fabry-Perot cavities, are represented by Dirac delta 
functions Q to constrain the k integration 

2 ^ 

^^o \l + il^l^)e^i-^^2r,\2 =Yd^^(k~k,,)' (13) 



where fc„ = TT[n+ {1 — 6<j)/TT)/2]/d, and rj —> m the limit 
of infinite potential strength. Here Scj) is the difference in 
the spinor polarizations of the two scattering potentials, 
and S(j) — 27m denotes the situation for identical scatter- 
ers. An incoming wave vector satisfying the resonance 
condition in Eq. gets fully transmitted through the 
two-barrier system. The modes in between the barriers, 
on the other hand, are fully refiected yielding the appro- 
priate quantization condition. Combining these results 
we obtain 



F = 2 



— k 
27r 



1 - J 5] '5(fc - 

n=0 



The Casimir force in Eq. (|14p can be calculated by 
applying the generalized Abel-Plana formula, 



tdt 



E 

n=0 



tdt 



sinh(27ri) 



Vcos(27r/3) — cosh(27rt) 



1 K 



(15) 



which is valid for < (3 < 1. Due to the rapid con- 
vergence of the integral in Eq. (fT51) , the result does not 



require an introduction of an explicit ultraviolet cutoff 
function More generally, since the refiection coeffi- 
cient vanishes at high energy it will regularize the calcu- 
lation of the force. Using Eq. (jl5p we obtain the force 
for two barriers satisfying bag boundary conditions. 



24d2 



1 



(16) 



for — TT < (50 < TT beyond which it is periodic. We also 
explore the force between two scatterers of finite height 
and width. In the small barrier strength limit the force 
becomes 



F2 cos{S(P) 



27rd2 

The force in the limits of F - 
plotted for three periods in 
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oo and F <c; 1 for a 
b in Fig. H 



(17) 
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FIG. 2: Force between two barriers as a function of their 
relative spinor polarization (50. The solid and dashed lines 
represent the forces in Eq. (|16|l and Eq. (I17p . respectively. The 
magnitude of the force in the F ^ 1 limit, the dashed curve, 
is rescaled to F = 1/2 so the two curves can be compared. 

The scaling of the force with distance as 1 /d'^ and the 
ratio of 1/2 between the repulsive and attractive forces 
are universal results for massless one-dimensional fluctu- 
ating fields in the limit d ^ a. When the range of the 
potentials becomes comparable to their separation, the 
first order correction due to the shape of the scatterer 
scales with 6F/F ~ a/d as seen in Eq. (|17p. analogous 
to a multipole expansion of an electrostatic interaction. 

The relative orientation can be expressed as S(j) — 
cos-i(I4 • l72/(|yi| . 11/2 1))- When the two potentials are 
aligned at (5(/) = 27rn we have F = — 7r/24c?2. This yields 
the attractive fermionic Casimir force as found in Ref. 

When Scj) = (2n -I- l)7r the relative polarization of 
the defect potentials is antiparallel and F = 7r/12c?2, i.e. 
a repulsive Casimir force is obtained. An analog of our 
result for a one-dimensional bosonic field is obtained by 
imposing mixed Dirichlet and Neumann boundary con- 
ditions where attractive and repulsive Casimir forces are 
found for like and unlike boundary conditions, respec- 
tively . A Casimir force that oscillates as a function of 
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defect separation d is known to arise from large momen- 
tum backscattering (Friedel oscillations) of the Fermi gas 
Q. However, the interaction we calculate here is mono- 
tonic as a function of distance. In our calculation, the 
magnitude and sign of the force varies as a function of 
the relative polarization of two scatters at a fixed dis- 
tance. As shown in Fig. [2] this behavior occurs for both 
finite barriers and hard- wall boundaries. 

The cusps seen in Fig. [2] at the odd multiples of n 
result from a sum over the discrete number of energy 
levels En {6(j)) . The energy bands found in Eq. (fT5)) cross 
zero energy at Scf) = {2n + l)7r as shown in Fig. [S] At 
fixed chemical potential, with negative energy states of 
the Dirac sea occupied, the number of states changes 
by one in each 27r periodic region indicated by dotted 
vertical lines in Fig. [31 Consequently, the force exhibits 
a discontinuity in slope in Fig. [2] exactly at the values of 
(j) at which there is a jump in the number of occupied 
energy levels. When the barrier strength is finite, the 
cusps in the force disappear. The resonance condition 
resulting in quantized states between the barriers is only 
valid for hard-wall boundaries. Note, the energy states 
between finite barriers exhibit a continuous spectrum. 
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FIG. 3: Quantized energy bands for massless Dirac fermions 
due to hard-wall boundary conditions as a function of the rel- 
ative polarization of the two potentials Sep. Solid lines denote 
the energy levels of the filled Dirac sea. Vertical dashed lines 
define 2tt periodic states where the number of occupied states 
changes by one. 

The interaction Eq. is likely to be important for 
defect interactions on carbon nanotubes, and possibly 
for other one-dimensional systems as well. Reinserting 
dimensional factors this force corresponds to an inter- 
action energy Ec = —Trhvp/'^^d for two identical scat- 
terers. With hvp ~ 5.4 eV- A this gives an energy of 
1.4 meV at a range d = 50 nm. Note that its spatial 
form follows the same scaling law as the Coulomb in- 
teraction between uncompensated charges, but it is re- 
duced by a factor 7r?it;_F/24e^ ~ .05. Thus, for charge 
neutral dipoles p — es whose electrostatic interactions 
scale as E^ ~ —p^/d^ = — (e^/d) x {s/dY, they are dom- 
inated by the Casimir interaction in the far field d > 5s. 



Similarly, this one-dimensional Casimir interaction com- 
pletely dominates the familiar van der Waals interactions 
between charge neutral species that are mediated by the 
fiuctuations of the exterior three dimensional electromag- 
netic fields. 

In order to fully understand the Casimir effect between 
defects on carbon nanotubes, one needs to consider the 
symmetry and range of the potentials produced by lo- 
calized defects. The spinor polarization discussed in this 
paper is determined by the form of the impurity poten- 
tial: fjz and Gy potentials define a sublattice-asymmetric 
and bond-centered defects, respectively. In addition, the 
electronic spectrum contains two distinct Fermi points 
at inequivalent corners of the two dimensional Brillouin 
zone. Short-range potentials couple the two Fermi points 
resulting in intervalley scattering [l3|. Therefore, both 
the structure of the defects and the effect of interval- 
ley scattering determine the sign and magnitude of the 
Casimir interaction. In the context of our model, a sharp 
potential is one with a range on the order of the tube ra- 
dius R for which the effects of intervalley scattering are 
suppressed by a factor of ad R, where ac is the width of 
the graphene primitive cell. Atomically sharp scatterers, 
on the other hand, will usually require a treatment of the 
effects of intervalley as well as intravalley scattering. 

To summarize, we introduced a force operator ap- 
proach for calculating the Casimir effect and obtained 
the fiuctuation-induced force between two finite square 
barriers mediated by massless Dirac fermions in one di- 
mension. In taking the limit of sharp barriers of infi- 
nite strength we obtained a Casimir force that scales as 
and is tunable from attractive to repulsive form as 
a function of the relative spinor polarizations of the two 
scattering potentials. 
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